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Abstract. Strange quark matter in beta equilibrium at high densities is studied in 
a quark confinement model. Two equations of state are dynamically generated for the 
same set of model parameters used to describe the nucleon: one corresponds to a chiral 
restored phase with almost massless quarks and the other to a chiral broken phase. 
The chiral symmetric phase saturates at around five times the nuclear matter density. 
Using the equation of state for this phase, compact bare quark stars are obtained with 
radii and masses in the ranges i? ^ 5 — 8 km and M ~ Mq. The energy per baryon 
number decreases very slowly from the center of the star to the periphery, remaining 
above the corresponding values for the iron or the nuclear matter, even at the edge. 
Our results point out that strange quark matter at very high densities may not be 
absolutely stable and the existence of an energy barrier between the two phases may 
prevent the compact quarks stars to decay to hybrid stars. 
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1. Introduction 

The recent discoveries of X-ray sources, whose origin may be attributed to strange 
stars is stimulating theoretical research on the structure, composition, dynamics 

and evolution of such objects. It is believed that the innermost region of compact 
stars, such as neutron stars, can be made out of hyperons, meson condensates or quark 
matter [H IHl 13 E] , whose thermodynamical behavior can be expressed by means of a 
relativistic equation of state (EOS) jHIHlini- 

Various effective models for the nucleon, using quarks as fundamental fields, have 
been used to provide EOS's for quark matter, which are subsequently applied to 
investigate the structure of compact stars [21 El HDl El- In this vein, we use the 
chromodielectric model (CDM) ^21 1^ to obtain EOS's for dense quark matter, and 
report, in this letter, the properties of the resulting quark stars. The chiral CDM 
provides a reasonable phenomenology for the nucleon [T^ IT^ and N-N* transition 
amplitudes ^HI- Baryons appear as solitons with three quarks dynamically confined 
by a scalar field, whose quanta can be assigned to 0"'^''" glueballs. In the quark matter 
sector, the model yields a relatively soft EOS at large densities [13 CHI- Drago et al. 
applied the CDM with a quadratic potential to the structure of compact stars [TTH OH]. 
They connected the EOS for quark matter, as provided by the CDM, with an EOS for 
hadron matter, obtaining stars with masses in the range 1 — 2Mq and radii of the order 
10 km with a hadron crust of 2 km. 

In the present work we consider an extension of the model used in Ref. ^Hj, now 
taking quartic instead of quadratic potentials. For such potentials two qualitatively 
different EOS's are dynamically generated for the same set of model parameters (the 
same used to describe the nucleon). They correspond to two distinct phases: in one 
phase the x fi^ld is large and the quarks are massless (chiral symmetric phase). In 
the other one, x is small and quarks are massive (chiral broken phase). The EOS's 
are degenerate at very high densities but with an energy barrier between them. As 
the density decreases, they separate and saturate at very different densities, yielding 
different energies per baryon number. The EOS for the chiral symmetric matter 
saturates at densities much higher than the nuclear matter equilibrium density and 
with an energy per baryon number above that quantity in iron or in nuclear matter at 
low densities. The compact stars obtained using this EOS are made out of quarks only, 
since the density at the edge is much above the nuclear matter saturation density, and 
hadronization processes do not take place. 

The prediction of smaller and denser objects in comparison with the neutron stars 
is quite exciting in view of the recent discovery of X-ray sources, by the Hubble and 
Chandra telescopes, which increased the plausibility that these sources might be strange 
quark stars [IIEIIS]- In particular, the isolated compact object RX J1856. 5-3574, not 
showing evidence for spectral lines or edge features pT| l22j. reinforced the conjecture for 
the existence of strange matter stars. However, the first results indicating a small radius 
for that compact object were based on the reported parallax of Ref. which was not 
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taking into account the camera geometrical distortion. Considering these corrections 
the parallax was reduced by roughly a factor of two [23 and the distance has 
been revised from 60 to 117 pc approximately. The revised radius {R = 11.4 ± 2 km) 
and mass [M/Mq = 1.7 ± 0.4) make that isolated compact object reproduceable by 
many EOS with and without strange matter [211 EH] • Hence, it is not guaranteed that 
RX J1856. 5-3574 might only be described as a quark star. 

2. Quark matter 

The CDM Lagrangian can be written as [T^ IT^ IT^ 

^ ~l~ ^a,n ~l~ -^g— meson ~l~ ' ("^) 

where 

jC.,. = ^d^ad^'a + Id^^ ■ d^^^ - Wiji, a) , (2) 

(vr, a) being the usual Mexican hat potential for the chiral mesons sigma and pi. The 
quark-meson interaction, assuming two flavours, is 

The last term in contains the kinetic and the potential piece for the x-field: 

C^ = \d,xd'x-U{x)- (4) 
We consider a quartic potential 

where M is the x mass. It has a global minimum at x = and a local one at x = • 
The height of the local minimum, B = [/(7M) = {rjM)'^, may be interpreted as a "bag 
pressure" [2Zj, as in the MIT bag model, and that will be used to fix parameters in 

u{x). 

In the soliton sector of the model, a good description of the nucleon is obtained for 
X close to zero, a region of the potential © where the cubic and quartic terms play no 
role. All dependences go into just one parameter, namely G = \/gM, and best nucleon 
properties are obtained for G ~ 0.2 GeV (e.g. M = 1.7 GeV and g = 0.023 GeV 
In non-strange homogeneous quark matter the deppendence on the single parameter G 
is even exact for the quadratic model [2HI- In our study we keep that value for G. In the 
quark matter sector of the model with double minimum potential there are solutions 
with small x ^.nd solutions with large i-e. in the region of the local minimum of the 
potential, as discussed below in more detail. From B = rj^M^ and assuming the wide 
range 0.150 < B^/^ < 0.250 GeV, one has 0.08 < r/ < 0.15, using M = 1.7 GeV. The 
parameter 7, which does not affect the phenomenology of the homogeneous matter, is 
not a totally free parameter: 7^ > 677'^ since the quartic term in ^ must be positive. 
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The extension of the model to include the strange quark, requires that one more 
term be added to the interaction Lagrangian Q, accounting for the coupling between 
the strange quark and the x field. We take it in the simplest form 

■^s— meson '4^s'^s ■ (6) 

X 

In addition, for the sake of beta equilibrium, an electron gas must also be included. 
The mean-field energy per unit volume for a homogeneous system of u, d and s quarks, 
interacting with x ^"^^ plus electrons, is 

d^k f— — : r^-^ d^k 







,2 



+ u{x) + ^,{<y'-fl)\ (7) 

[f-K = 93 MeV and we always use m^r = 1.2 GeV in this paper, though the results are 
not sensitive to the sigma mass]. The first two terms refer to quarks, and the third 
one to the electrons, all described by plane waves. The degeneracy factor is a = 6 (for 
spin and color). The Fermi momentum for each type of particle, ki, is related to the 
corresponding density, pi, through pi = akf/^Gn"^). The quark masses in are all 
different [23 : 

gu,d (7 Qs 

Xfn X 

with the coupling constants for each flavour given by Qu = g {fn + ^3), gd = g {f-K — ^3) 
and gs = g — fn) [(,3 = —0.75 MeV, fx = 113 MeV]. Since the vacuum expectation 
value of the conflning fleld is zero, the quark masses raise up to inflnity for densities 
approaching zero. 

A variational principle applied to the energy density, Eq. ((Tj), leads to two gap 
equations for a and x- In the interior of a compact star the matter should satisfy both 
the electrical charge neutrality and chemical equilibrium. These conditions together 
with the gap equations, form a system of six algebraic equations solved at each baryon 
density p = | {Pu + Pd + Ps)- For a given p one obtains a self-consistent set a, x, k^, kd, 
kg, ke and, with such solution for each density, the energy density [Eq. 0] is readily 
evaluated as well as the energy per baryon number, e/p, the pressure, etc. 

For the same set of parameters we found two distinct stable solutions, hereafter 
denoted by I and II. For both solutions, a remains always close to /^r irrespective of the 
density. In solution I, the x held, close to zero, is a slowly increasing function of the 
density. For small the quartic potential ^ is indistinguishable from U = |M^x^, 
thus, in practice, our solution I is the one obtained and used by Drago et al. in 
the framework of the quadratic potential. Due to the smallness of the x held, quark 
masses are large [see Eq. (jHJ] and the system is in a chiral broken phase. On the 
other hand, for solution II, the conflning fleld is large, x ~ 7^ (local minimum of U), 
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almost independent of the density. The resulting quark masses are similar for the three 
flavours and very close to zero (chiral restored phase). Because the chemical potentials 
are dominated by the Fermi momentum, one has fiu — f^d — f^s, and therefore fie — 0, 
i.e. in solution II there are almost no electrons. Besides solutions I and II, there is an 
additional unstable solution corresponding to x ~ 7M/2 [local maximum of U{x)]- 

For each solution we obtained the corresponding energy per baryon number as a 
function of the baryon density (EOS) (see Fig. Q). EOS-I is not sensitive to 7 and rj 
(since x is small), just depends on G, and it is rather similar to the one used in Ref. |19j . 
The saturation density occurs at a low density, slightly higher than the nuclear matter 
equilibrium density, po- Its shape, at intermediate densities, is similar to hadronic EOS's 
(see Ref. ^7] for the two flavours sector). 




Figure 1. Energy per baryon number versus density for solutions I (solid line, small 
x) and II (dashed line, x ^ jM) and various parameters rj (other model parameters: 
M = 1.7 GeV, g = 0.023 GeV and 7 = 0.2). The dotted hue corresponds to the 
unstable solution with x ~ 7M/2. 

The EOS-II is also insensitive to 7, but does depend on rj [in fact, the dependence is 
on {rjM)'^, as we have already discussed]: the energy per baryon number increases with 
rj and so does the saturation density. Depending on rj the minimal energy per baryon 
number of solution II can be either below or above solution I. For 77 ~ 0.12 the saturation 
occurs at p ~ 5po and the energy per baryon number is some 230 MeV higher than for 
solution I at its saturation density. The two stable solutions are almost degenerated 
at high densities in the narrow range 0.1 < rj < 0.12. In Fig. ^the dotted lines refer 
to the EOS for the unstabe solution of the gap equations (supplemented by electric 
neutrality and beta equilibrium conditions), for which the x is at the local maximum 
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of the potential (0)- In order to undergo a transition from I to II, the system has to 
go through the energy barrier represented by the dotted EOS. The barrier gets higher 
at small densities and, for rj = 0.12 (third panel in FiglH), at p ~ 5po, e/p ~ 2.7 GeV 
for the unstable solution. Therefore a transition from one regime to the other is not 
likelly to occur and both minima in the EOS I and II are stable. In a 3D plot of the 
energy per baryon number versus (p, x) the stable solutions correspond to two distinct 
"valleys" , and the unstable solution corresponds to the top of the barrier between these 
two valleys as it is shown in Fig. El Our results point out that the metastability of 
strange quark matter only occurs at high densities for solution II. 




Figure 2. 3D plot of the energy per baryon number versus density and x field. 
Solutions I lies in the valley for small x, the unstable solution is at the maximum and 
solution II lies in the valey at high x- The parameters arc as in Fig. ^ 

Interesting enough, our results with two distinct EOS, are consistent with the results 
from a recent calculation in perturbative QCD [20]. We checked that, for 77 ~ 0.12 the 
CDM reproduces accurately the EOS's obtained in Ref. jHOl- However, in our case, there 
is no parameter fit to get two equations of state: they are dynamically generated for the 
same Lagrangian parameters. Also in Ref. jHI], using an extension to finite chemical 
potential of lattice QCD data for the equation of state, similar results were obtained. 
Therefore, there seems to be some model independence, which is worth to point out 
here. 

To investigate the structure of stars we solved the Tolman-Oppenheimer-Volkoff 
equation. Since EOS-I is identical to the one using a quadratic potential, it leads to 
stars that have the same phenomenology as the hybrid stars obtained by Drago et 
al. pnj: -R ~ 10 — 12 km, a hadron crust and a mass M ~ 1 — 2Mq. 
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Since EOS-II saturates at a high density and, in addition, the system is not hkely 
to undergo a transition to solution I, one should not perform any connection to the 
hadronic sector: the EOS-II alone generates a new family of strange quark stars. In 
Fig.Elit is shown the mass-radius relation for different values of 77. These quark stars are 
smaller and denser in comparison with those resulting from EOS-I. For rj ~ 0.115 (and 
M = 1.7 GeV, yielding B^^^ ~ 0.195 GeV) one obtains a maximum radius i? ~ 6 km 
(and a corresponding mass M ~ O.QM©). According to our calculation, such star has a 
central density of lOpo (po is the nuclear matter density) and a central energy density 
e ~ 3x10^^ g/cm^. At the edge, the density drops to 5po and e ~ 1.35x10^^ g/cm^. The 
ratio e/p remains approximately constant inside the star and the minimum period of 
the star, computed using the expression given in Ref. [22]; is ~ 0.4 ms. 




Figure 3. Mass versus radius for the pure quark stars (solutfon II) in the CDM model. 
For the other model parameters see caption of Fig. 

A possible mechanism to explain the formation of such bare strange stars could be 
a supercooling effect in the early universe, when a significant amount of quark matter 
was frozen into quark stars. This mechanism needs an effective potential with two 
minima, as discussed in [33] , resembling the quartic x potential. As we see from Fig. El 
the mass-radius relation for these strange small stars mainly depends on the height of 
the local minimum of the x potential. Finally, it is also worth mentioning that our 
results for the maximum radius and maximum mass for the stars from solution II are in 
agreement with Ref. jHl], where it was shown that the Chandrasekhar limit for quarks 
stars depends on B. 

Let us summarize our results. We have considered the CDM with a quartic 
potential, with parameters fixed in the nucleon sector and to yield a reasonable 
bag constant. Using a mean- field variational method we obtained two solutions for 
homogeneous strange matter in beta equilibrium, one similar to the already known 
solution for quadratic potentials (with massive quarks) and a new one with massless 
quarks. The pure quark stars emerging from the chiral symmetric solution are small 
and dense compact objects. To get so a small stars, the saturation density should be 
high enough, and this is achieved with a corresponding energy per baryon number always 
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above the energy per nucleon in the iron at low density or the energy per nucleon in 
nuclear matter at its saturation density. 
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